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On the existence of hyperplane sequences, 
with quality parameter and discrepancy bounds 

Friedrich Pillichshammer and Gottlieb Pirsic * 

Abstract 



It is well-known that digital (t, m, s)-nets and (T, s)-sequences over a finite field 
have excellent properties when they are used as underlying nodes in quasi-Monte 
Carlo integration rules. One very general sub-class of digital nets are hyperplane 
nets which can be viewed as a generalization of cyclic nets and of polynomial lattice 
point sets. In this paper we introduce infinite versions of hyperplane nets and call 
these sequences hyperplane sequences. Our construction is based on the recent 
duality theory for digital sequences according to Dick and Niederreiter. We then 
analyze the equidistribution properties of hyperplane sequences in terms of the 
quality function T and the star discrepancy. 



1 Introduction 

In the field of quasi-Monte Carlo (QMC) integration one approximates the integral of a 
function over the (usually high-dimensional) unit-cube by the average of function eval- 
uations over a well-chosen deterministic point set. Such an integration rule is called a 
QMC-algorithm. The fundamental error estimate for this method is the Koksma-Hlawka 
inequality which states that the absolute error of a QMC-algorithm is bounded by the 
variation in the sense of Hardy and Krause of the integrand times the star discrepancy 
^ 1 of the underlying point set. Hence, point sets from the unit-cube having a very uniform 

distribution, or in other words, having low star discrepancy yield small integration errors, 
at least for functions with bounded variation in the sense of Hardy and Krause. For more 
information on QMC we refer to the books [31 18|. 

The star discrepancy of a point set & = {x , . . . ,x^_x} consisting of N points in 
[0, l) s is defined as 

D* N (&) = sup - X S (J) 

where the supremum is extended over all intervals of the form 

J= [0,a x ) x ... x [0,a s ) C [0,l) s , 

where A(J,N) is the number of indices n for which x n belongs to J and A S (J) is the 
s- dimensional volume of J. For an infinite sequence S = {y n ) n >o f° r N G N the star dis- 
crepancy Dtf(S) is the star discrepancy of the point set consisting of the first iV elements 
of S. 
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Currently the best constructions of point sets with low star discrepancy are based on 
the concept of (t, m, s)-nets in base b as introduced by Niederreiter [7J (see also [3] and 
[HI Chapter 4] for surveys of this theory). 

Definition 1 ((£, m, s)-nets). Let b > 2, s > 1 and < t < m be integers. A point set 
& consisting of b m points in [0, l) s forms a (t,m, s)-net in base b, if every subinterval of 
the form rLLiM~ d % ( a * + l)b~ di ) of [0, l) s , with inte gers di > and integers < < b 1 
for % = 1, . . . , s and of volume 6*~ m , contains exactly b l points of 

In the definition above, t is often called the quality parameter of the net. It is clear 
from the definition that any point set consisting of b m elements from [0, l) s is at least an 
(m, m, s)-net in base b. Smaller values of t imply stronger equidistribution properties for 
nets. 

Infinite versions of (t, m, s)-nets are (t, s)-sequences as introduced by Niederreiter [7J 
and their generalizations (T, s)-sequences as introduced by Larcher and Niederreiter [6]. 
We present here a revised version introduced in [15] and occasionally called 'digital se- 
quences in the broad sense' that eliminates certain problems with ambiguous digit expan- 
sions. For this we introduce the truncation operator 



Note that the definition is dependent on the given expansion. We also extend this operator 
to act on vectors by component-wise application. 

Definition 2 ((T, s)-sequences). Let b > 2, s > 1 be integers and let T : No — > No be a 
function which satisfies T(m) < m for all m G No- A sequence S = (xq, x±, . . .) of points 
in [0, l) s forms a (T, s)-sequence in base b, if for all m, k e No, the point set consisting of 
the points [xp.b m ]b,m, ■ ■ ■ , [%(k+i)b m -i}b,m forms a (T(m), m, s)-net in base b. 

In the definition above, T is often called the quality function of the sequence. Any 
sequence in [0, l) s is at least a (T, s)-sequence in base b with T(m) = m for all m G No- 
Smaller values of T imply stronger equidistribution properties for sequences. A (T, s)- 
sequence in base b is called a strict (T, s)-sequence in base b if the value of T(m) can not 
be decreased at least by 1 for any m G No- It is known from [3l Theorem 4.32] that a 
strict (T, s)-sequence in base b is uniformly distributed modulo one if 



If there exists a t G No such that T(m) = t for all m G No then one speaks also of 
(t, s)-sequences. Obviously, (t, s)-sequences are uniformly distributed modulo one. 

Explicit constructions of (t, m, s)-nets and (T, s)-sequences are usually based on a 
digital method over a finite field which will be explained in the following: 

From now on let q be a prime-power and let ¥ q be the finite field of order q. Let 
<p : {0, 1, . . . , q — 1} — > F q be a fixed bijection with ip(0) = 0. 

Definition 3 (digital (t, m, s)-nets over ¥ q ). Let s,m G N. Let C\,...,C S be m x m 
matrices over ¥ q . For i = 1, . . . , s and for k G {0,l,...,g m — 1} with g-adic expansion 



in 




i=l 



i>l 



lim m — T(m) = oo. 
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k = k + Kiq + • • • + K m -\q m 1 with Kq, ■ ■ ■ , K m -l G {0, 1, . . . , q — 1} multiply the matrix 
Ci by the vector k = (<p(kq),. . . , ^(/? m _i)) T , i.e., 

Qk=: (y tA (k),...,y hm (k)) T G (F™) T , 

and set 

._ p^iviM) , <p~ 1 (yi, m (k)) 

If for some integer t with < t < m the point set consisting of the points 

x k = (xk,i, ■ ■ ■ ,x k ,s) for k G {0, 1, . . . ,q m - 1} 

is a (i, m, s)-net in base q, then it is called a digital (t,m, s)-net over ¥ q , or, in brief, a 
digital net. The matrices Ci,...,C s are called its generating matrices and the matrix 
C = (CjCj . . . Cj) G F^ xsm is called the overall generating matrix of the digital net. 

The construction of infinite sequences follows more or less the same lines as above, 
more detailed: 

Definition 4 (digital (T, s)-sequences over F 9 ). Let s G N. Let C lf . . . , C s be N x N 
matrices over ¥ q . For i = 1, . . . , s and for k G No with g-adic expansion k = kq + 
Kiq + ■ ■ ■ with k ,Ki, . . . G {0, 1, . . . , q — 1} multiply the matrix Ci by the vector k = 
O(/t ),y?(Ki),---) T > i.e., 



and set 



Gk=: (y M (A;),^ 2 (A;),...) T G(Fj)~ 



q q A 



If for some function T : No — > No with T(m) < m for all m G No the sequence consisting 
of the points 

x k = (%k,i, Xk,s) for k G N 

is a (T, s)-sequence in base q, then it is called a digital (T, s)- sequence over ¥ q , or, in 
brief, a digital sequence. The matrices Ci, . . . , C s are called its generating matrices. 

It is known from [21 Theorem 4.86] that a strict digital (T, s)-sequence in base b is 
uniformly distributed modulo one, if and only if 

lim m — T(m) = oo. 

Many constructions of digital nets are inspired by a close connection between coding 
theory and the theory of digital nets (see, for example, Niederreiter [91 [101 [11] ) . Examples 
for that are the so-called (u, u + v) -construction (see [U [13]), the matrix-product con- 
struction (see [12]) and the Kronecker-product construction (see [U E]). Here we deal 
with a construction for digital nets which is an analog to a special type of codes, namely 
to cyclic codes which are well known in coding theory. This construction has been intro- 
duced by Niederreiter in [9] who adopted the viewpoint that cyclic codes can be defined by 
prescribing roots of polynomials. Later this construction has been generalized by Pirsic, 
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Dick and Pillichshammer [18] to so-called hyperplane nets whose definition will be given 
in Section [3] where also some results on the quality parameter are recalled. 

It is the aim of this paper to introduce a corresponding construction also for digital 
sequences and to prove some results concerning the equidistribution of the resulting se- 
quences. Such a construction can be obtaind with the help of duality theory for digital 
sequences as recently introduced by Dick and Niederreiter [2]. 

The paper is organized as follows: In Section [2] we recall the basic facts from dual- 
ity theory for digital nets and sequences. This theory will be an important tool in our 
construction. Furthermore, in Section [3] we recall the definition of and some basic results 
on hyperplane nets. The construction of what we call hyperplane sequences will be pre- 
sented in Section HI In this section we also show how to determine the quality function of 
hyperplane sequences and we prove an 'if and only if condition under which hyperplane 
sequences are uniformly distributed modulo one. In Section [5] we determine the generator 
matrices of hyperplane sequences and in Section E] we compare hyperplane sequences to 
Larcher- Niederreiter sequences. Last but not least, in Section [TJ we show the existence of 
hyperplane sequences which satisfy T(m) = slog^m + 21og g logm + 0(1). This in turn 
implies the existence of hyperplane sequences which satisfy a certain bound on the star 
discrepancy. 



Duality theory is an important tool for the construction of digital nets and sequences and 
provides a connection to the theory of linear codes. It has been introduced by Niederreiter 
and Pirsic [H] for digital nets and by Dick and Niederreiter (2] for digital sequences. We 
also refer to [3, Chapter 7] for an overview. In the following we recall the basic facts which 
will be necessary for the construction of hyperplane sequences in Section HJ 



2 Duality Theory 



2.1 Duality Theory for Digital Nets 

Let be a F g -linear subspace of ¥ s q m . Then its dual space AT C F q m is its orthogonal 
complement relative to the standard inner product in F^ m , i.e., 



{A e 



<sm 



B 



A = for all B 6 Af}. 



Definition 5 (NRT weight). For a 



(ai, . 



. , a m ) e F™ let 




if a = 0, 
aj ^ 0} if a ^ 0. 



We extend this definition to F^ m by writing A e ¥ s q m as the concatenation of s vectors of 
length m, i.e., A = (a 1; . . . , a s ) E F* m with 6 F™ for i = 1, . . . , s, and putting 



s 




i=l 



The weight V m is called the Niederreiter-Rosenbloom-Tsfasman weight. 
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Definition 6 (minimum distance). Let M ^ {0} be a F^-linear subspace of ¥ s q m . Then 
the minimum distance of Af is defined by 

5 m (N) = min{y m (A) : A 6 jV\ {0}}. 

Furthermore, we put 5 m ({0}) = sm + 1 by convention. 

The following theorem has been first shown by Niederreiter and Pirsic [T3] . It can be 
also found as [3], Corollary 7.12]. 

Theorem 1. Let m, s G N, s > 2. Then, from any ¥ q -linear subspace M of ¥ s q m with 
dim(A/") > sm — m one can construct generating matrices Ci,...,C s of a strict digital 
(t, m, s)-net over ¥ q such that 

t = m- 6 m (AT) + 1. 
The row space C of the overall generating matrix C satisfies C = Af ± . 

2.2 Duality Theory for Digital Sequences 

A digital sequence over ¥ g is fully determined by its generating matrices C 1: . . . , C s G 
F^ xN . For m G N we denote the m x m left-upper sub-matrix of Ci by C\ . The 

matrices C[ m \ . . . , are then the generating matrices of a digital (t, m, s)-net over ¥ q . 
The overall generating matrix of this digital net is defined by 

for any m G N. 

Hence a digital sequence can equivalently be described via the sequence 
of overall generating matrices or the sequence Ci, C2, . . . of row spaces thereof. Each 
or C m C Ff" describes the first q m points of the digital sequence and each has a dual 
space C ¥ s q m (which is the null space of in F^ m ) associated with it (see [31 

Section 7.1]). Hence the dual for a digital sequence consists of a sequence (C^) m >i of dual 
spaces which have certain relations to each other, that we shall explain in the following, 
crucial definition. 

Definition 7 (dual space chain). Let s G N, s > 2. For all m G N, let M m be a F^-linear 
subspace of F* m with codim(A/" m ) < m. Let 

A/" m+ i, m := {(a lr .,a s )6F; ra : ((a 1; 0), . . . , (a s , 0)) G Af m+1 } , 

i.e., J\f m +i im is the projection of vectors in Af m+ i with zeroes at every (m+l)-th coordinate. 
Suppose that Af m+ i tm is a F g -linear subspace of Af m with dim.(M m / Af m +i,m) — 1 f° r an 
m G N. Then the sequence (J\f m ) m >i of spaces is called a dual space chain. 

The following result has been first shown by Dick and Niederreiter [2]. It can also be 
found as [31 Corollary 7.27]. 

Theorem 2. For a given dual space chain (Af m )m>i, one can construct generating ma- 
trices Ci, . . . , C s of a strict digital (T, s) -sequence over ¥ q such that 

T(m) = m - 8 m {H m ) + 1 for all m G N. 

For all m G N the mth row space C m of the mth overall generating matrix satisfies 
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3 Hyperplane Nets 



Hyperplane nets are an example of digital nets whose construction is based on duality 
theory. They are a generalization of cyclic nets as introduced by Niederreiter [§] and were 
first presented in [18]. The slightly more general definition given here was first presented 
in |17j . An overview of hyperplane nets can be found in [3] Chapter 11]. 

Definition 8 (hyperplane nets). Let m, s £ N, s > 2, and a prime-power q be given. Let 
/ £ ¥ q [x] with deg(/) = m and let R m := ¥ q [x]/(f) be the polynomial residue class ring 
modulo the ideal (/). For i = 1, . . . , s let Q5j = {bj,i, . . . , bi, m } be an ordered basis of the 
polynomial residue class ring R m considered as vector space over ¥ q . 

Fix an element a. = («i, . . . , a s ) £ R s m \ {0} and consider the subspace 

W a := {k £ R s m : a ■ k = (mod /)}. 

Define the mapping 8 : R s m — > ¥ s q m by 

k = . . . , k s ) £ R s m i-> . . . , /ti, m , . . • , /c a ,i, . . • , K a , TO ) £ F* m , 

where . . . , «j,m) is the coordinate vector of fcj £ i? m with respect to the basis QSj. 
We denote by Cq, the orthogonal subspace in F^ m of the image M a := d{N' a ). Let 

C a £ F™ xsm 

be a matrix whose row space is C a . Then we call the digital net with overall generating 
matrix C a a hyperplane net over ¥ q with respect to Q3i, . . . , *B S . This hyperplane net will 
be denoted by & a and we say £P a is the hyperplane net associated with ex. 

Note that the dual space of a hyperplane net is 

N<x = 0(K) = {A £ ¥ s q m : C a A = 0}. 

In the following we show how to determine the quality-parameter t of a hyperplane 
net. To this end we need a further definition: 

Definition 9 (figure of merit). For a £ R s m the figure of merit is defined as 

s 

p{a) = s - 1 + min V deg(^) = -1 + 5 m {M a ). (1) 

kGA^\{0} ^ 

The following result gives a relation between the quality-parameter of a hyperplane 
net and its figure of merit. It can be found as [16, Theorem 1]. 

Theorem 3. A hyperplane net & a associated to a £ R s m \ {0} is a digital (t,m, s)-net 
over ¥ q with t = m — p{at). 

Observe that there is no dependence on the chosen bases. This means that a different 
choice of bases can be considered as a scrambling of the original sequence that does not 
deteriorate the quality. 

The following existence result is a slight generalization of [HI Theorem 2]. 
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Theorem 4. Let m, s G N, s > 2 and /et q be a prime-power. Choose ordered bases 
Bi, . . . , B s of R m over ¥ q . For p G Z define 



S-l ✓ s P+Ci , 

d=0 ^ ' 7=0 ^ 



/ s — d + 7 — 1 

7 



<f + 1 - q p+s 



g-i 

9 



s-l 



, i/jen t/iere exisi more t/ian 



Let ft G R snc/i that < ft < 1. IfA q (s,p) < /3q m (s 

(1 — ft)((q — l)g m_1 ) s_1 elements a. = (1, a2, . . . , a s ) £ -Rm gcd(aj,a;) = 1 for all 
i = 2, . . . , s, such that p(oc) > s + p. The associated hyperplane nets are digital 
(t, m, s)-nets over ¥ q with t < m — s — p. 

Corollary 1. Let m, s G N, s > 2 and m sufficiently large, and let q be a prime-power. 
Let ft G R such that < ft < 1. There exist more than (1 — j3)((q — l)g m_1 ) s_1 elements 
a = (1, «2, • • • , a s ) G R s m with gcd(«j, x) = 1 for all i = 2, . . . , s, such that 



P{<*) > 
Proof. For p > 1 we have 



(s - I V 

m + log 9 ft-(s- l)(log, m - 1) + log 9 — 



^ <- g(3(.-^(:!::i)S 



d=0 

< <f+'£ 



d=0 

s-l 



_ ir - d -i q Ap+s-iy- d - 1 



(s-l) 



s-d-l)\ 



where O s indicates that the implied constant depends only on s. Now let 

(«-!)! x 



m + \og q ft -(s-l) \og q m + \og q ( 



which is in N for sufficiently large m. Then 



A q (s,p) < fa* 



q-1 



1 + 



log q ft 



m 



-(s-l) 



\og a m 



+ — los 



J 9 e— 1 

m m q s 1 



s-l 



x ( 1 + O s ( - 



g-l 



s-l 



for sufficiently large m. The result follows from Theorem HI 



□ 



4 Definition of Hyperplane Sequences 

For m G N set f m (x) = x m G F g [x]. For each i = 1, . . . , s and all m G N choose a basis 
25m,i — {^m,i,i, • • • , &m,i,m} of the polynomial residue class ring R m := ¥ q [x]/(x m ) over ¥ q 
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such that, expressed in the canonical basis, {1, x, x 2 , . . . }, we have b m +i,i,j = (bm,i,j, 0) for 
all 1 < j < m. Define the mapping 9 m : R s m — > W q m by 

k = [ki, . . . , k s ) G R m I—)- (ACl,l, • • • , Kl,mi ■ ■ ■ i K s,li ■ ■ ■ i K s,m) £ F q j 

where (/c^i, . . . , «i, m ) is the coordinate vector of fcj G i? m with respect to the chosen basis 
<8 mj ; for % = 1, . . . ,s. 

Let F g [[x]] be the space of formal power series over ¥ q . Let Y denote the subset of 
power series whose reductions modulo any x m are invertible modulo x m , i.e., 

Y = {a G ¥ q [[x\] : gcd(a mod x m , x m ) = 1 Vm G N}. 

Let a e 7 s and, for m G N, set (mod x m ) (coordinate- wise). 

Theorem 5. For a G Y s and m G N let W {m) = {k G R s m : a (m) • k = (mod x m )}. 
Then the sequence (J\f a 

( m ) )ra>l ; 

where U a ( m ) = 6 m (M' a(m) ), is a dual space chain. 

Proof. Let Mm+i,m be the set 

Af m+ i, m ■= {(a!,..,a s )eF™ s : ((ai, 0), . . . , (a s , 0)) G A/" a(m+1) } . 

By the choice of the sequence of bases (23 m ,j) m >i for i — 1, . . . , s, it follows that jV m +i, m 
is a Fg-linear subspace of M a (m) . It remains to show that 

d\m(U a (m) / N m +i,m) < 1 dim(A^ m+ i, m ) > dim(A^ ( m)) - 1. 

We have 

M rn+1 , m = 9 m ({k G R s m : «( m+1 ) • k = (mod x m+l )}). 

Now we have A G M m+ i, m if and only if A = m (k) where c*( m+1 ) • k ee (mod x m+1 ). 
Let a( m+1 ) = (ai, . . . , a s ) G -R^+i where a, = a^o + cti,i% + • • • + a i:m x m for i = 1, . . . , s 
and let k — (ki . . . , fc s ) G i?^ with fcj = /c ij0 + k^x + • • • + Ki^ m _\X m ~ x for i = 1, . . . , s. 
Then a( m+1 ) • k ee (mod x m+1 ) is equivalent to 



E 

3=0 



\ 



EE 

= 1 l,u>0 
1+11=3 



V 



x 3 = 0. 



(2) 



/ 



Hence M m +i m is isomorphic to the nullspace of the (m + 1) x sm matrix A over F 9 given 

by 



E 

i=i 

Therefore we obtain 

dim(jV m+ i )m ) > sm — m — 1 = dim(jV a (m) ) — 1 
and the result follows. □ 



/ 


«i,0 







. 


\ 






«i,0 





. 






Oti,m—l 


<^i,m-2 


Cj,m— 3 


• CXi,0 




\ 




&i,m—l 






/ 
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Combining Theorem [2J, Theorem [5] and §Q we obtain the following corollary. 

Corollary 2. From any a. EY S one can construct a strict digital (T, s) -sequence S a over 
¥ q where 

T(m) = m - p(a (m) ) for all m G N. 

Definition 10 (hyperplane sequence). The strict digital (T, s)-sequence S a over F 9 is 
called a hyperplane sequence associated with a. 

Theorem 6. Let a = (ai, . . . ,a s ) G Y s . The hyperplane sequence S a over¥ q is uni- 
formly distributed modulo one, if and only if ax, ■ ■ ■ ,oc s are linearly independent over 
V q [x]. 

Proof. By [31 Theorem 4.86] a strict digital (T, s)-sequence over ¥ q is uniformly dis- 
tributed modulo one, if and only if linim^oom — T(m) = oo. Hence, by Corollary [2J the 
sequence S a is uniformly distributed modulo one, if and only if 

lim p(atS m ^) = oo. 

Let S a be uniformly distributed modulo one. Assume that there exist p%,...,p s G 
F 9 [x], not all of them the zero polynomial, such that 

pit*! H \-p s a s = G Fj[x]]. 

Then for m > maxi<j< s deg(pj) we have (pi, . . . ,p s ) G N\ m ) \ {0} and hence 

s 

p( a M)< s -l + ^deg(p i ). 

i=i 

Thus, p(o;( m )) is bounded as m tends to infinity which contradicts the uniform distribution 
of S a . 

Let a. = (ax, . . . , a s ) G Y s be such that ai, . . . , a s are linearly independent over ¥ q [x]. 
Assume that the sequence S a is not uniformly distributed modulo one. Then there exists 
a K G N such that 

p(a (m) ) < K for all m G N. 

From the definition of the figure of merit p it is evident that a bound on p implies a 
bound on the degrees of the possible pi for a nontrivial linear dependence relation. Thus 
it follows by the pigeonhole principle that there exist pi , . . . , p s G ¥ q [x], not all of them 
the zero polynomial, such that 

p ia [ m) + ■ ■ ■ + p s a^ ] = (mod x m ) 

for infinitely many m G N. This however leads to a contradiction, since by the linear 
independence of ax, ■ ■ ■ , oc s over ¥ q [x] it follows that pi«i + • • • +p s a s ^ G F 9 [[#]]. Hence 
there exist Z G No and G F 9 for i > I and a; 7^ 0, such that 

Pi«i H h p s a s = aix 1 + ai +1 x l+1 H 

and therefore 

Pxa^ + ■■■+ p s a { ™ ] £ (mod x m ) for all m > I. 

□ 
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5 Explicit Generator Matrices 



First we briefly recall the explicit form of hyperplane net generator matrices (cf. [IB] and 
[21 Section 11]). From this, the shape of the infinite generator matrices of hyperplane 
sequences will then be evident. 

Going back to the general ring R m = ¥ q [x]/ (f) of Definition [HJ an injective ring homo- 
morphism \l/ : R m — > F™ xm can be constructed by setting equal to the companion 
matrix of the generating polynomial f(x) and, for all further elements, requesting the ring 
homomorphism property, Wx,y,z G R m : *&(xy — z) = *ff(x)ty(y) — ^(z). The validity of 
this construction can be realized using the fact that the minimal polynomial of a com- 
panion matrix is the polynomial itself. Also, it is easily verified that this homomorphism 
condition is equivalent to setting 



' m— 1 



m— 1 



Va 1; ...,a m G Fq : * ^ 



x 



. i=0 



i=0 



Next, take the canonical basis (£j = {1, x, . . . , x m ~ 1 } for R m , as a F g -space. For i — 1, . . . , s 
let Bi be the matrix associated to the basis change from the given basis 53 j to (B{. Then 
the following was proven in [18] (the version given here follows [3j Theorem 11.5]; the 
theorem as originally stated in [18] had a minor flaw): 



Theorem 7. The matrices Ci := (^(ct^-Bj) 
matrices for the hyperplane net & a . 



s can be chosen as generating 



associated to the 
for 



Now note that in the case of hyperplane sequences the matrices B, 
bases 23 m ,i are upper triangular: this follows inductively, since the > 
m > and i = 1, . . . , s consist of the coordinates of bm,ij, j = 1, . . . ,m with respect to 
the canonical basis and we requested the last component to be equal to unless j = m. 
Let 7r m denote the action of taking the upper left mxm submatrix of a finite or infinite 
matrix. Then it also follows inductively that for m' < m we have ^m'{B mi ) = B m > ti . 
Hence there are infinite upper triangular matrices Bi G F^ xN such that 7r m (i?j) = B m i . 

Denote with \]>( m ) the ring homomorphism as discussed above, associated to the ex- 
tension ¥ q [x]/(x m ) used in the construction of the hyperplane net & a ( m ). It can be given 
explicitly quite easily: first note that the companion matrix of the polynomial x m is the 
matrix with entries of 1 in the first lower subdiagonal and elsewhere, so ^^(x) = 



(4 

consequently 



i J+1 )i j'=0,...,m— 1 • 



It follows that ^ m \x k ) = (5 iij+k ) 



i,j+k)i,j=0,...,in- 



_i for k = 0, 



m 



1 and 



( a ii0 
a i,l a-ifi 

\0'i,m— 1 (^i,m—2 ■ ■ ■ 



■ o\ 

. 



where a 



(m) 



E 

j>0 



SO 



\]>( m )( 



a 



(m)x 



is a lower triangular Toeplitz matrix. We see that each of these matrices is 



the upper left mxm submatrix of an (also lower triangular Toeplitz) infinite matrix in 



pNxN -yyhjQh we denote as ^/(a,;) 



Thus we have shown that the generator matrices of <^ Q (m) can be chosen as C m ^ : = 
(7r m (\I/(aj))7r m (-Bj)) T . However, by the triangular structure of the matrices also the prod- 
uct ^(a)Bi is admissible and we have furthermore C mi = 7r m (\l/(aj)_Bj) T . This implies 
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that the first block [cco]&,ro> ■ ■ ■ , [ x b m -i]b, m of the sequence & a coincides with £P a ( m ). We 
summarize: 



Then, the generator matrices of & a are given by Ci = (^(a^Bi) 1 . 

Remark 1. In particular cases, e.g., if the Bi have finite columns (or, equivalently, all bij 
are polynomials) all generator matrices have the finite-column property, i.e., almost all 
entries are 0. This is the property required for the definition of 'digital sequences in the 
narrow sense'. Hence, if the bases include nonpolynomial Laurent series, and consequently 
the matrices would lack this requirement of finite columns, then for sequences in the 
'narrow' sense we would only be able to state: for m > 0, the generator matrices of 
7T m (^ Q ( m )) can be chosen as C m ^ = 7r m (\l'(Q; i )_Bj) T . 

Remark 2. Note that Ci = B T ^(ai) T is the product of a lower triangular matrix with an 
upper triangular Toeplitz matrix and can therefore be interpreted as a linear scrambling 
(in the sense of Faure and Tezuka [I]) of the matrices \l/(aj) T . 



6 Relation to Larcher-Niederreiter Sequences 

Recall that Larcher-Niederreiter sequences, as defined in [HI Ch. 4.4, Eq. 4.64] and treated 
in detail in (5] are based on elements chosen from ¥ q ((z)) (where z := 1/x), or rather, 
actually the subring zFjfz]] is used for the definition. It is a natural question to ask what 
the relation to hyperplane sequences is. 

The shortest way to describe Larcher-Niederreiter sequences is by defining them as 
being generated by arbitrary infinite Hankel matrices; in the construction cited above, 
the entries of a matrix map bijectively to the positive-degree coefficients of a Laurent 
series. It is well-known that finite Hankel matrices correspond to polynomial lattice point 
sets (PLPSs), see [21 Chapter 10]. Hence initial 6-adic blocks of Larcher-Niederreiter 
sequences are, up to truncation, equivalent to PLPSs. 

It has been shown in [T7J Theorem 2] that PLPSs are also a special case of hyperplane 
nets, when choosing the canonical basis: This choice results in the Bi being identity 
matrices. So in the case of the <^ a (m) (with the canonical basis) for any m > we have 
Cm t i = \E , (o; J - m ' ) ) T as generator matrices. Let p( m ) be the horizontally reflected identity 
matrix. Then the matrices C mi iP( m > still generate the same net, since the regular matrix 
p( m ) only effects a reordering of the point set, but the matrices are now (triangular) 
Hankel matrices, hence generate a PLPS. 



Theorem 8. Let a £ Y s and let *& m ,i be as in Section^ For i = 1, . . . ,s define the 
matrices Bi £ ]F^xN ^ setting, for j £ N, the first j entries of the jth column equal to 
the coefficients of bj ti j modulo x J and zero elsewhere. For ai = Yl,j>o a i,j x ^ e ^MM] f or 
i — 1, . . . , s associate infinite lower triangular Toeplitz matrices ^(oti) £ xN by 




Qi,m— 1 ^i,m— 2 • • • ®>i,0 
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We will show in the following that, however, the matrices C$ = "lim m C m ,i" do not 
generate a Larcher-Niederreiter sequence (to motivate this intuitively, consider that, as 
m grows, the matrices p( m ) shift the first few points of <^ a (m) ever further to the end of 
the sequence, effectively removing them in the limit). For our proof we first restate [21 
Theorem 3] (cf. also [31 Th. 7.28]) and add a corollary. 

Theorem 9. Let (Af m ) m >i be a dual space chain. The generating matrices Ci, ...,C S are 
unique up to a multiplication ofCi, C s from the right with the same nonsingular upper 
triangular (NUT) matrix, if and only if dim(Af m +i,m) = sm — m — 1 for all m G N 

Corollary 3. If a sequence is generated by regular matrices Ci, . . . ,C S that fulfill 

rank(C (m ' m+1) ) = rank(C (m) ) + 1 = m + 1 for all m G N, (3) 

where C^ m+1 ^ := (cj m ' m+1)T , . . . , ci m ' m+1)T ) and the cj m ' m+1) are the left upper m x 
(m + 1) submatrices of Ci, then any s matrices generating the same sequence have the 
form C\ = CiP, with P an arbitrary NUT matrix. 

Proof. This follows from Theorem H] by noting, as in Section [21 that M m is the dual space 
of C m , the row space of C^ m,m \ Due to the regularity C m has dimension m, hence its 
dual M m has dimension sm — m. Furthermore, J\f m +i, m is the dual space of the row 
space of C^ m ' m+1 \ Hence, if the rank condition is fulfilled, the dimension of A/" m +i. m is 
sm — rank(C ( '" v ' n+1 )) = sm — m — 1 and the theorem can be applied. □ 

Now what is left is to verify the condition fl3]) for the C m ,i as given above. Then, by the 
corollary, any matrices generating the same sequence would have to be upper triangular, 
hence are never Hankel matrices as would be the case for Larcher-Niederreiter sequences. 

Note that the first m row vectors of £^ m > m+1 ) are ^ ne m _ dimensional coefficient vectors 

of x^af 1 ' (mod x m ) and the final, (m + l)-th row vector is equal to (a- m+1 ^ — ai^)/x. 
Supposing a nontrivial linear combination of the row vectors of (7( m . m + 1 ) exists, then, for 
some Xj G ¥ g , j = 0, . . . , m — 1 and for each i = 1, . . . , s we get 

m— 1 



XjX^oq 71 = (a J - m+1 ' ) — a^oj/x (mod x m ) 

3=0 

m— 1 

<^ XjX 1 a i = (a i — ai t o) fx (mod x m ) 

3=0 
m—l 

<^ AjX" ,+1 a i = ctj — a i)0 (mod x m+1 ) 

3=0 

^A (m) (x)a, = ai,o (mod x m+1 ), 



for some nonzero polynomial A^ m \x). But this can only occur if all a ™ are constant 
multiples of each other. Requesting the «j to be linearly independent over ¥ q , which is 
the only interesting case by Theorem [61 rules this out for large enough m, i.e., at least 
for large m the condition (J3} of the corollary is satisfied. But this suffices, since the 
proof of Theorem actually shows how to construct C^ m+1 ^ uniquely, given 
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Furthermore, any other C G xN generating the same sequence would have to fulfill 
c ,(m) = c M p (m) ) • = wit h pH being a NUT matrix. But this can never 

happen, if the Cj are also NUT matrices and the C[ are Hankel matrices. Hence Larcher- 
Niederreiter sequences and hyperplane sequences have at most trivial common generator 
matrices and are distinct as constructions. 



7 An Existence Result 

In this section we prove an existence result which is based on Corollary [TJ In the following 
let 

Y s := {a = (a!, . . . ,a s ) <EY S : a x = 1}. 

Theorem 10. Let (3 = ((3 m )m>i be such that ^ m>1 m < oo. Then there exists an integer 

L = L(s, q, (3) and an ct £Y S such that the associated hyperplane sequence S a is a digital 
(T, s)- sequence over¥ q with 

s — 1 

m 

T(m) < log 9 — h L for all m > 1. 

Pm 

For example if we choose (3 such that f3 m = m( - lci g m ^ for m > 2, then we obtain the 
following corollary: 

Corollary 4. There exists an integer L and an at e Y s such that the associated hyperplane 
sequence S a is a digital (T, s) -sequence over ¥ q with 

T(m) < s log q m + 2 log 9 log m + L for all m > 2. 

Corollary 5. There exists a a G Y s such that the star discrepancy of the associated 
hyperplane sequence S a satisfies 

n*t<!\ ^ (logAQ 2s (loglogiVn f 77Ar> o 
D N {S a ) = 01 — J for all N > 3, 

where the implied constant depends only on s and q. 

Proof. Combine Corollary H] with the star discrepancy bound given in [3] Eq. (5.5), p. 
196]. □ 

For the proof of Theorem [10] we need some notation: Let Y m = {a (m) : a G Y}. Note 
that \Y m \ = q m ^ 1 {q — 1). We define now a probability measure over the set Y . We do 
this in a way such that the measure of corresponding vectors in Y m is equiprobable. For 
m G N let /i m be the equiprobable measure on the set Y m . We say a subset A of Y is of 
finite type, if there exists an integer m = m(A) G N and a subset A' of Y m such that 

A = {a G Y : a {m) G A'}. 

The measure of the finite type subset A is then defined as 

\A'\ 

fi(A) = fl m (A') = qm _ 1{q _ 1 y 

13 



(Of course, the number m = m(A) is not uniquely defined by A. However, it is easy to 
see that the definition of /i does not depend on the specific choice of m.) 

For s G N let /jl^ denote the complete product measure on Y s induced by fi. 

Now we are ready to prove Theorem [TU1 

Proof. Let 

U m (P) = {a G Y s : p( a ^) > m + log, - (s - l)(log, m - 1) + log, |. 

Then it follows from Corollary CD that there exists some mo = mo(s,q) such that for all 
< f3 < 1 and all m > m we have 

n.(u m (P))>i-p. 

Let U m {0) be the complement of U m {(3). Let {3 = (/3 m ) m >i be such that ^ m>1 /3 m < 
oo. Let n = n (s, q, (3) be such that n > m and J2 m > no ^ m — ^- Then we have 



f| W m (/3 m ) J = l-/z a I |J W ro (^ m ) 

Vm>no / \m>rao 

> 1 - J2 VsiUmiPm)) 



m>no 

> i - E « 

> 



TO 

m>no 



Hence there exists aaGT such that 
p(« (m) ) > 



0-1 

m + log P m -(a- l)(log m - 1) + log 



<? g s-l 



for all m > n . 



From this it follows that there exists some integer L = L(s, q, f3) and there exists a ex. G Y s 
such that 

p(a^) > |_m + log, (3 m — (s — 1) log, m — Lj for all m > 1. 
Now the result follows from Corollary [2J □ 
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